Collective versus single-particle effects in the optical spectra of finite electronic 

quantum systems 
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We study optical spectra of finite electronic quantum systems at frequencies smaller than the 
plasma frequency using a quasi-classical approach. This approach includes collective effects and 
enables us to analyse how the nature of the (single-particle) electron dynamics influences the op- 
tical spectra in finite electronic quantum systems. We derive an analytical expression for the low- 
frequency absorption coefficient of electro-magnetic radiation in a finite quantum system with bal- 
listic electron dynamics and specular reflection at the boundaries: a two-dimensional electron gas 
confined to a strip of width a (the approach can be applied to systems of any shape and electron 
dynamics - diffusive or ballistic, regular or irregular motion). By comparing with results of nu- 
merical computations using the random-phase approximation we show that our analytical approach 
provides a qualitative and quantitative understanding of the optical spectrum. 



Optical spectra of finite metallic systems have been 
intensively investigated for almost a century. Early ap- 
proaches such as Mie's Q are of classical nature. In |l]] 
the absorption of electro-magnetic radiation by conduct- 
ing spheres is determined. It is shown that the absorp- 
tion spectrum exhibits a resonance at tu p /y3 (where u) p 
is the bulk plasma frequency) due to collective motion of 
the charge carriers. 

In the last decades there has been a substantial amount 
of work on the nature of such collective resonances in 
metallic clusters 0, nuclei ||, thin films small 
metal particles [p] , and dimcnsionally reduced quan- 
tum systems fp|-[l0f investigating, in particular, quantum- 
mechanical effects. In most of these cases, the elec- 
tron dynamics is ballistic (the mean free path £ is larger 
than the system size a). The majority of theoretical pa- 
pers on the Mie resonance of finite metallic systems use 
the so-called random phase approximation (RPA), a self- 
consistent, quantum-mechanical approach incorporating 
collective effects. The nature of the Mie resonance in fi- 
nite electronic quantum systems is well understood, both 
qualitatively and quantitatively (for a recent review see 
for instance ||). 

The emergence of the field of Mesoscopic Physics has 
fueled an increased interest in the electronic and op- 
tical properties of finite, disordered quantum systems 
(with diffusive electron dynamics, f < o) in external 
fields. In this context, attention has largely focussed on 
(quasi-)static properties, to ~ A/h, where A is the mean 
level spacing of the system in question (which is gen- 
erally smaller than tu p by many orders of magnitude). 
Mesoscopic fluctuations of the static polarisability and 
the capacitance, for instance, were characterised in [ pr| , 
and the electron density itself in In all of these 

cases, collective effects must be taken into account in or- 
der to adequately deal with the screening of the external 
field. In the static limit, a Thomas- Fermi (TF) ansatz is 
appropriate [ pr|Jl^ ]. 



Much less is known about optical spectra in the low- 
frequency region A/h < uj <C w p , despite the fact that 
this regime is of particular interest: one expects that 
the spectra strongly depend on the nature of the (single- 
particle) electron dynamics: In ballistic systems, for ex- 
ample, it was argued E| that optical spectra should ex- 
hibit resonances near multiples of u> c — irv-p/a (see also 
ft4|l). Usually A/h < lu c < uj p . According to [jj], these 
resonances should overlap to give a linear frequency de- 
pendence of the absorption coefficient for ui S> w c in two- 
dimensional systems (classical and local electro-magnetic 
theories predict a quadratic dependence independent of 
dimension). These arc striking and unexpected results. 
They were, however, obtained within a TF approxima- 
tion which is valid in the static limit. It must be ex- 
amined to which extent dynamic screening effects may 
modify the results Jl5| . 

Furthermore, in systems with ballistic chaotic dynam- 
ics, the optical spectra should, at least within a single- 
particle picture, reflect universal energy-level correlations 
found in classically chaotic quantum systems (ll|. This 
issue was addressed in the seminal paper by Gorkov and 
Eliashberg []l7|], investigating the polarisability of metal- 
lic particles with disordered walls. However, the treat- 
ment neglects screening effects all together jl8|-^0| . How 
does dynamical screening modify this picture? 

Numerical calculations (based, for instance on the 
RPA) are ill-suited to answer these questions: They pro- 
vide little qualitative insight in the low-frequency region 
and, more importantly, it is necessary to consider small 
systems or to make use of symmetries in order to make 
the numerical computations feasible. Disordered and 
(asymmetric) chaotic systems are very difficult to deal 
with. In order to understand how the classical single- 
particle dynamics is reflected in optical spectra of finite 
quantum systems, it is thus greatly desirable to have an 
analytical theory incorporating collective effects. 

The main result of this paper is an analytical expres- 
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sion for the absorption properties for a two-dimensional 
electron gas confined to a strip of width a [eqs. ( |l4|) and 
(|l5|)]. This example exhibits all features of more compli- 
cated geometries. However, in the case discussed here, 
numerical RPA calculations are feasible and allow us to 
discuss the accuracy of the analytical approach. 

In the following it is assumed that A/h < uj <C 
Lu p ,Ep/h where E-p is the Fermi energy. It is further- 
more assumed that A, A s a 3> Ap where A is the wave 
length of the external radiation, A s is the skin depth and 
Ap is the Fermi wave-length. 

We consider a closed metallic quantum system placed 
in an external electric field -E ox t [with a time dependence 
tx exp(icut)]. If the wave-length is much larger than the 
system size a, the external field is approximately constant 
throughout the system and (neglecting retardation ef- 
fects) can be written as the gradient of an electric poten- 
tial (/?ext( r ) — E a xe x [r = (x, y, z) is a three-dimensional 
coordinate vector]. Within the RPA, the electronic re- 
sponse of the system to E cxt is calculated by solving a 
set of self-consistent equations for the effective electrical 
potential ip{r) = (p cx t(r) + 5ip(r) [Sf(r) is the potential 
due to the induced charge density Sg(r)] 



Stp(r)= J dr'G(r,r')6g(r'), (1) 
Sg(r) = — J dr'n (r, r 1 ; uj) <p(r) 

with the boundary condition that Sip(r) vanishes as 
\r\ — > oo. G(r,r') is the Green function of the Laplace 
equation AG(r,r') = — e~ a 1 5(r — r'), e is the dielectric 
constant. Ho(r,r';uj) is the non-local polarisability 



n (r,r'; W ) = -2e 2 ^ 



a,/3 



/(fa) ~ /(gff) 

E a — ep — tbu+vy 



xKWMr'WpWMr) , (2) 

e is the electron charge, e a and ip a ( r ) are the single- 
particle eigenvalues and eigenfunctions of the undis- 
turbed system, they are usually calculated with in a 
Hartree-Fock or a local-density approximation, /(e) = 
&(Ep — e), and 7 > is smaller than A. Within the 
RPA, the absorption coefficient (proportional to the en- 
ergy dissipation per unit time) may be written as H 



huj 



where 



d(u) 



J drdr'Sg* (r)H (r,r';uj)(p(r') 



(3) 



(4) 



is the (complex) dipole moment and the asterisk denotes 
complex conjugation. 

In the following we derive an explicit analytical ex- 
pression for the absorption coefficient a(u>) of a finite 



electronic quantum system in an external electric field, 
valid in the frequency range A/h < uj C Wp. We make 
use of the approximations suggested in BE. First, ac- 
cording to Fermi's golden rule, the absorption coefficient 
is 



2 A 2 El 



(ip a \(p\ipp) 



±E F 
■£0 r - 



(5) 



Second, the matrix elements of <p are evaluated within a 
semi-classical approximation JT^ , p2| -p5|. Third, <p itself 
is determined within a quasi-classical approximation: ac- 
cording to (Q), the effective electric potential is given by 
(in symbolic notation) 



(p = -II 1 5g. 



(6) 



In order to determine Sg, eqs. (|l|) are usually solved nu- 
merically, using a real-space discretisation [^6| or by ex- 
panding in a suitable basis set. An approximate analyti- 
cal solution may be obtained by noting that for uj <C oj p , 
||GIIo|| S> 1. In other words, 8g is well approximated by 
the classical charge density Sg c i of the metallic system 
subjected to an external potential ip e xt, 



Atp c i = -Sg c i/e 



(7) 



with tp c i(r) — > ip cx t{f) as \r\ — > 00 and ip c \ = within 
the system (in the classical limit, the external field is thus 
screened out completely) . Eq. (0) may be solved for 8g c \ 
using standard methods |27[]. Fourth, LTo is determined 
within a quasi-classical approximation [[28|]2l| 



n (r,r';cj) = e 2 v d 5{r - r') + iujP {d \r,r'- uj) (8) 

where v& is the density of states per unit volume in d 
dimensions and P^ d '{r,r';ui) is the Fourier transform of 
the classical propagator (r,r';t). In ballistic systems 
it is written as a sum over classical paths p from r to r' 



pW(r,r» =2 

cl. paths p 



det 



d(r') 



d(T p ,n p ) 



exp(kjTp) . (9) 



Here t p is the time taken from r to r' along the path 
p, and rip is a unit vector describing the direction of the 
initial velocity. For diffusive systems see e.g. [f29f . 

In the following we show by comparison with quantum- 
mechanical RPA calculations that provide a qual- 
itative and quantitative description of absorption in the 
frequency range A/h < u -C u) p . 

Two-dimensional strip. We consider a two-dimensional 
electron gas confined to a strip in the a:-y-plane sur- 
rounded by vacuum J30|, subject to a time-dependent 
electric field E cxt directed along the negative x-axis [com- 
pare fig. |l|(a)]. The width of the strip (along the x-axis) 
is a, its length L (along the y-axis), with L 3> a. Within 
the system, the electrons move ballistically, and they are 
specularly reflected at the boundaries at x = ±a/2. 
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We write Sg(r) = 5<r(x,y) 5(z). For L 3> a, the 
surface-charge density 8a depends on x only and the re- 
solvent G is written as 



e Q J 2tt 2\q\ 



(10) 



With 



- J dydy' n (r; r'; u>) = A (x, x'; lo) S(z) S(z') , (11) 

the RPA-equations (|l|) are reduced to a set of one- 
dimensional equations for ip(x, z — 0) with the kernel 
Ao(x, x'\ lo). We model the confinement in the x-direction 
by introducing hard-wall boundary conditions. This is 
adequate in the range of parameters considered below 
and simplifies the quasi-classical analysis. We solve the 
resulting self-consistent equations numerically using a 
real-space discretisation and obtain the absorption co- 
efficient from (||). 

The corresponding quasi-classical approximation for 
a(u>) is obtained as described above: the classical 
surface-charge density is determined by solving (M) in 
elliptic cylinder coordinates: 5a c \(x) — 2eoEo(a 2 /'i — 
x 2 )~ 1/>2 <d (\x\ — a/2). The corresponding classical field 
lines are shown in fig. |l|(a). According to eqs. (|8|J^) 
the one-dimensional kernel Ao(x, x' '; lo) is given by a sum 
over classical paths as shown in fig. [y(b) which may be 
summed by Poisson summation. Using (^|) one obtains 
for the effective electric potential 



ip(x) — cos [fin (x/a + 

ju>0 



1/2)] 



(12) 



e 2 v 2 



(?Wf) 



Eq. (|T^) has an intuitive interpretation: neglecting a 
correction term (phdy(%) which is small except for x in 
a boundary layer of width 8x oc vf/lo, eq. (12|) can be 
written as a sum over two terms 



V — Vstat + Vdyn 



(13) 



where Lp s ta,t{x, z — 0) — (eo q s )~ 1 5cr(x) is the (linearised) 
TF potential (q s — e 2 v 2 jeo is the two-dimensional TF 
screening vector) and fdyn is a dynamical contribution 
corresponding to a current building up the screening 
charges. It obeys d 2 ipd yn /dx 2 = — (m e Lo 2 /e 2 ) Sa/ao 
where m e is the electron mass and Co is the areal charge 
density of the electrons. 

Fig. ||(a) shows tp(x) according to ( p"2| ) and ( [l3| ) com- 
pared with the results of a numerical RPA calculation. 
One observes excellent agreement (and </?bdy is small 
except at the boundary). Our results show that for 
larger frequencies (u> > u> c ), the dynamical potential 
fdyn makes a significant contribution to ip and dynamical 
screening effects cannot be neglected. For the absorption 
coefficient we obtain using eqs. (||) and (jlj) 



. . TT 2 he?,a w 2 



e 2 v 2 



odd 

E 

H>[u/u a ] 



VP 



-J X 2 W2) (14) 



with limiting forms 

, \ _ j ^Che^auJ 2 / {e 2 U2io c ) for w<oj c , 
a \ UJ ) — <y 7r fi e ^aui/(4:e 2 V2) for lu > lu c 



(15) 



and C ~ 0.12. In fig. g(b), we show quantum-mechanical 
RPA results in comparison with eqs. (|lj) and (15) and 
observe excellent agreement. We observe prominent res- 
onances in the absorption coefficient near odd multiples 
of lo c , due to single-particle cyclotron orbits (electrons 
moving in phase with the external field). This estab- 
lishes that the single-particle resonances conjectured in 
|l3| exist within the RPA. Their positions, strengths and 
shapes are very well described by (jlj). Eq. (|l5| ) and the 
inset of fig. @(b) show that in the limit of large frequen- 
cies (to ^> lo c ), the absorption coefficient is linear in lo. 
In the opposite limit (lo <C lo c ) where the TF approach 
is adequate, it is quadratic. 

We conclude that the quasi-classical approximation de- 
scribed above, for the parameters considered here, pro- 
vides a quantitative description of the optical properties. 

Three-dimensional thin film. To conclude we discuss 
ip for a thin film [|TJ of width a in the y-z-plane sub- 
ject to an external potential tp ext (r) = Eoxe x . The 
classical charge density is concentrated at the boundary, 
5g c \(x) — ±<5er §{x ± a/2). The corresponding static and 
dynamical contributions to ip{x) are singular; we thus use 
the TF charge density || instead of 8g c \ (appropriate in 
the limit of small k s corresponding to high electron den- 
sities): SgTF{x) = k s eoE smh(k s x)/ cosh(k s a/2). Here 
k 2 = e 2 u^/eo is the three-dimensional TF screening vec- 
tor. The RPA equations are easily solved within a real- 
space discretisation. Our results for ip(x) are shown in 
fig. ||, and compared to results of the analytical approach 
using eqs. (||]). 

We have also calculated the absorption coeffcient for 
A/h < lo <C Lo p within the RPA. The analytical approach 
must be used with caution in the case of the film since 
it requires that ip be smooth on the scale of Af- It turns 
out that (x(lo) is, to a good approximation, quadratic in 
lo as opposed to the two-dimensional case. As in d = 2 
dimensions we observe resonances near odd multiples of 
lo c (not shown). 
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FIG. 1. (a) Electric field lines (in the x-z-plane) for an 
infinitely long metallic strip of width a (in the :r-y-plane, ori- 
ented along the y-axis) placed in a constant external electric 
field .Eext = —Eoe x . (b) Classical paths from x to x' con- 
tributing to Ao(x, x'; lj). 




FIG. 2. Left: <p(x) for a strip of width a — 10 [a.u.] with 
r B — 1: quantum-mechanical results (•), analytical results ac- 
cording to eq. ^ ( ) and eq. ^ ( ). The inset 

shows the correction term ipbdy(x). Right: Shows Imd(w) 
(for a = 10 4 [a.u.] and r B = 1) as a function of ui: RPA result 
( ) and eq. (114) ( ). The inset shows a(cj) as a func- 
tion of lo: RPA result ( ) and eq. ( |lE| ) ( ). As usual 

r s = ro/ao where ao is the Bohr radius and ro is the length 
scale defined in terms of area per electron. 




FIG. 3. Shows (p(x) for a thin film of width a = 50 [a.u.] 
with r s = 1: RPA results (•), analytical results according to 
eq- @ ( ) and using tp ~ ip 8tat + ip dyn ( ). 
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